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FOREWORD 


The  program  described  here  consisted  of  determination  of 
turbulent  atmospheric  heating  rates  and  diffusion  coefficients  through 
microdensitometric  reduction  and  computer  analysis  of  photographs  of 
chemical  releases.  Related  earlier  work  by  the  PhotoMetrics  research 
group  is  described  in  References  1,  2,  and  4. 

The  author  wishes  to  express  his  thanks  to  Mr.  Samuel  P. 
Zimmerman  (Technical  Monitor)  of  AFGL  for  his  continued  encourage- 
ment and  support,  to  Dr.  P.  I.  Richards  for  his  assistance  in  modify- 
ing AFGL  atmospheric  codes,  and  to  Dr.  I.  L.  Kofsky  and  C.  C.  Rice 
who  also  made  important  contributions  to  the  studies  described  in 
this  report. 
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SECTION  I 


INTRODUCTION 

The  objective  of  the  program  described  here  is  the  determination 
of  atmospheric  heating  rates  and  turbulent  diffusion  coefficients  from 
transport  and  surface  brightness  distribution  of  tracer  chemical  clouds, 
and  incorporation  of  these  values  into  AFGL  model  atmosphere  codes. 
Radiance  distributions  which  give  mass  densities  in  these  artificial  clouds 
are  recorded  photographically,  accessed  by  microdensitometric  scanning, 
and  analyzed  by  digital  computer  methods.  The  primary  approach  to 
this  problem  has  been  that  of  analyzing  the  "turbulent"  component  of 
spatial  brightness  distributions  of  releases  after  removing  low  fre- 
quency trends  by  filtering.  The  method  employs  Fourier  analysis  to 
determine  wavenumber  spectra  of  the  higher-frequency  fluctuations  in 
scene  brightness.  These  spectral  features  that  characterize  turbulence 
are  then  analyzed  to  determine  heating  rates  (rate  of  dissipation  of 
turbulent  kinetic  energy)  and  diffusion  coefficients  as  described  in 
Section  II.  Additional  estimates  of  heating  and  diffusion  have  been 
derived  from  measurements  of  horizontal  wind  fields  in  the  strato- 
sphere and  mesosphere  obtained  by  photographic  triangulation  methods 
(Section  III).  Work  with  simple  AFGL  model  atmosphere  codes  is 
described  in  Section  IV. 

A summary  of  chemical  release  cloud  data  analyzed  is  shown 
in  Table  I.  The  results  of  these  and  other  analyses  were  reported 
as  they  were  completed  to  AFGL  scientists  for  their  immediate  use 
in  interpreting  atmospheric  turbulence  phenomena. 


Figure  1. 


Turbulence  present  in  photograph  of  release  Joan 
(Wallops  I.  - 1974)  from  approximately  90-105  km. 
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SECTION  II 


DETERMINING  ENERGY  SPECTRA,  HEATING 
RATES,  AND  EDDY  DIFFUSION 

The  methods  by  which  turbulent  energy  spectra  may  be  calculated 
have  been  treated  in  Reference  1.  Briefly,  it  is  possible  to  charac-  1 

terize  turbulent  flow  in  the  atmosphere  (Ref's  2-5),whether  shear 
dominated  or  buoyancy  dominated,  by  analyzing  the  surface  brightness 
distributions  of  sunlight-scattering  or  chemiluminescent  tracers 
released  into  the  atmosphere  from  rockets.  After  a short  period  of 
time  the  chemical  reaches  pressure  equilibrium  with  ambient  species 
and  conforms  to  the  flow  field.  Photometrically  calibrated  photographs 
of  the  release  freeze  the  spatial  characteristics  of  the  flow  in  time  and 
according  to  the  ergodic  hypothesis  a scan  across  this  frozen  flow 
field  should  be  statistically  equivalent  to  the  temporal  output  of  a fixed 
physical  sensor  as  the  flow  proceeds  across  it.  A representative 
photograph  of  such  turbulence  near  100  km  (Release  Joan  - Wallops  I., 

1974)  is  shown  in  Fig  1. 

The  actual  measurements  are  made  by  microdensitometrically 
scanning  calibrated  photographic  images  of  the  glow.  These  photo- 
graphic densities  may  then  be  converted  to  scene  brightness  through 
the  H&D  characteristic  of  the  film  and  various  corrections  may  be 
applied  to  account  for  exposure  duration,  lens  T-stop,  and  varying 
sky  background  (Ref  2). 

Precise  photogrammetric  calibration  must  also  be  accomplished 
(see  next  section)  to  determine  the  altitude  at  which  the  measurements 
are  performed  and  the  range  to  the  cloud  must  be  known  to  convert 
film  plane  measurements  to  distances  at  the  glow.  For  analyses 
requiring  time  sequences  of  photographs,  correction  must  be  made 
for  range  variation  due  to  cloud  movement  by  the  mean  wind  field. 

Turbulent  energy  spectra  are  then  calculated  from  the  corrected 
brightness  data  either  by  standard  Fourier  transform  methods  or  by  a 


simple  summation  process  ("counting"  method  - Ref  7)  which  deter- 
mines the  energy  in  wavenumber  bands.  This  method,  which  is  also 
described  in  Ref  1,  calculates  the  energy  as  the  sum  and  average  of 
the  mean  squares  of  all  data  fluctuations  of  length  t multiplied  by  f . 

The  length  f is  defined  as  the  distance  between  zero  crossings  of  the 
high  pass  filtered  signal  and  estimates  of  power  may  be  calculated  at 
wavenumbers  k = 7r/f  . This  technique  is  very  rapid  and  provides 
reliable  spectra  in  the  mid  to  high  wavenumber  ranges  of  interest. 
Problems  with  data  truncation  and  filtering  end  effects  are  minimized 
(Ref  8).  It  does  require  a large  amount  of  input  data  but  this  is  also 
necessary  with  standard  Fourier  analysis  in  order  to  achieve  a desired 
statistical  accuracy  (Ref  9). 

Whichever  method  of  calculation  is  used,  correction  is  performed 
for  the  transfer  function  of  the  microdensitometer  aperture  and  for  the 
modulation  transfer  function  of  the  microdensitometer  itself,  which 
is  dependent  upon  temporal  scanning  rate  of  the  instrument.  It  should 
also  be  noted  that  true  microdensitometer  optics  (strict  control  of 
light  in  the  transmission  microphotometer,  and  spatial  stationarity  of 
its  response)  is  needed  to  ensure  sufficiently  accurate  photographic 
photometry  for  these  analyses.  Flare  (scattered)  light  and  nonuniform 
sensitivity  over  the  field,  present  in  most  full -field-illuminated  or 
"scanner"  type  optical  systems,  may  attenuate  or  nonuniformly  reproduce 
some  of  the  spatial  frequencies  present  (Ref  10).  Examples  of 
turbulence  spectra  calculated  for  release  Joan  are  shown  in  Figures 
2 and  3. 

Determination  of  heating  rates  from  the  spectral  (and  statistical) 
computations  may  now  be  accomplished.  Specifically,  the  atmosphere 
frequently  shows  a so-called  Kolmogoroff,  or  inertial  spectral  range  - 
a slope  indicating  minus  5/3-power  dependence  - that  implies,  from 
dimensional  arguments,  that  turbulent  energy  is  being  fed  at  a constant 
rate  from  large  cells  into  small  ones  in  the  subrange  of  cell  dimensions 
in  which  this  slope  obtains.  Whenever  this  inertial  subrange  is 
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present  the  following  arguments  may  be  utilized.  The  scalar  rate  of 
decay  (en)  of  the  density  fluctuations  (n')  due  to  molecular  diffusion  at 
a given  altitude  is 


where 


,n-6D<  ig)Z  > 

z r°° 

<(f^')  * J k2  E(k)  dk 


(1) 


and  D = molecular  diffusion  coefficient  dn'/dx  = space  rate  of  change 
of  turbulent  density  fluctuations  of  the  tracer  material  in  the  release, 
k is  the  wavenumber,  and  E(k)  is  the  energy  spectrum.  dn'/dx  may 
also  be  determined  from  the  brightness  data  by  numerical  differentia- 
tion. The  rate  of  viscous  dissipation  of  turbulent  kinetic  energy  is 


e 


A f 


E3  (k^  kt5 


(2) 


where  wavenumber  k^  is  located  within  the  inertial  range  of  the 

spectrum  and  A is  a constant.  This  method  of  determining  e has  been 

used  successfully  at  altitudes  near  100  km  (Ref  5 and  11),  but  is  limited 

to  ~80  km  and  higher  by  a combination  of  limited  photographic  spatial 

resolution  and  the  altitude  variation  of  the  kinematic  viscosity  v.  Theinte- 

2 1/3 

grand  in  the  equation  2 approximation  for  < (dn’/dx)  >,  behaves  as  a k 
to  ~0.  2 k^  for  the  inertial  range  of  E(k),  where  k^  = (e  /v  ) 4 is  the 
limiting  wavenumber  for  viscous  dissipation.  Thus  if  optical  resolu- 
tion limits  the  spectrum  E(k)  to  values  less  than  0.  2 k the  estimates 

2 v 

of  <(9n'/3x)  >,  and  c will  be  low.  In  existing  photographs  this  will 

begin  to  occur  roughly  at  80  km  and  grow  progressively  more  severe 

as  v decreases  with  altitude.  The  use  of  statistically  determined 

9n' /9x  generally  underestimates  these  quantities  as  well  since  a portion 

of  the  data  representative  of  sky  background  is  retained  in  the  calcu- 

2 

lation.  This  relatively  smooth  data  will  bias  the  average  of  (dn'/dx) 
toward  lower  values. 


/ 


12 


■SSB 


Another  approach  to  calculating  e may  be  taken,  which  goes  as 
follows.  The  average  of  the  square  of  the  turbulent  density  fluctuations 
n'  is  calculated  for  several  times.  Then 


€ 

n 


3(P'2) 


at 


(4) 


and  the  rate  of  turbulent  dissipation  is  calculated  from  equation  (3). 
This  formulation  for  decaying  turbulence  should  be  applicable  to  tracer 
pulse  or  trails  in  any  altitude  region. 

Determinations  of  eddy  diffusion  may  be  accomplished  utilizing 
the  turbulent  heating  rates.  It  has  been  shown  (Ref  6)  than  an  eddy 
diffusion  coeficient  may  be  expressed  as 


K = Cc 


(5) 


where  C is  a constant  and  kQ  is  the  low  wavenumber  limit  of  the 
"inertial"  range  of  the  turbulence  spectrum.  kQ  may  usually  be 
determined  by  inspection  of  the  turbulence  spectra  or  calculated  from 


k 
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oo 
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E (k)  dk 


r°°-i 

k E (k)  dk 
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or  estimated  as  Cq  (e  /N  ) for  buoyancy  dominated  turbulence 
where  is  the  length  scale,  Cq  = constant,  and  N = Brunt  Vaisala 
frequency. 


Further  measurements  of  heating  and  turbulent  diffusion  may  be 
made  from  knowledge  of  the  atmospheric  wind  field  (Section  III),  and 
diffusion  may  also  be  determined  by  using  configuration  space  methods 
similar  to  those  used  for  finding  molecular  diffusion  coefficients  by 
measuring  the  rate  of  growth  of  individual  cells.  Lastly,  low  frequency 
energy  spectra  may  be  used  to  calculate  wavenumber  spa  e halfwidths 
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of  clouds  which  retain  some  Gaussian  character  (recall  that  the 
Fourier  transform  of  a Gaussian  is  also  a Gaussian).  This  last  method 
may  however  show  very  large  errors  in  wavenumber  halfwidths  and 
derived  spatial  halfwidths  due  to  the  photographically  limited  data 
length  (ard  therefore  poor  resolution  of  low  frequencies). 


SECTION  III 


PHOTOGRAPHIC  TRIANGULATION  AND  VELOCITY  TERMINATION 

In  order  to  determine  quantitative  measures  of  turbulence  it  is 
necessary  to  accurately  know  the  spatial  location  of  the  chemical 
release  as  a function  of  time.  Over  the  time  of  a series  of  photo- 
graphic frames  required  to  calculate  diffusion  coefficients,  the  wind 
field  in  the  atmosphere  moves  the  cloud  with  respect  to  the  camera. 

This  change  of  range  introduces  a change  in  the  magnification  of  the 
image,  and  a corresponding  change  in  transverse  scale  or  conversion 
from  image  distance  to  dimensions  on  the  object  itself.  It  is  necessary 
to  triangulate  to  the  diffuse  glow  to  find  its  range  and  the  altitude  at 
which  the  diffusion  is  being  measured.  Methods  of  determining  these 
photo grammetric  calibration  factors  have  improved  in  accuracy  over 
a period  of  several  years  from  plane-trigonometric  to  spherical 
trigonometric  to  vector  solutions  (Ref's  12-16),  which  provide  altitude, 
latitude,  and  longitude  of  points  spaced  along  the  trail. 

The  triangulation  process  involves  determining  the  triangulation 
camera  orientation  from  photographs  of  the  stellar  background  as 
described  in  Ref  12,  Camera  position  in  latitude,  longitude,  and 
altitude  above  the  reference  ellipsoid  must  be  known  accurately  as  well 
as  the  time  of  the  exposure  which  records  the  star  positions.  Using 
this  time  and  the  known  camera  position  vector  equations  which  relate 
the  measured  film  plane  coordinates  of  stars  to  the  actual  star  positions 
in  right  ascension  and  declination  may  be  solved  to  determine  the 
orientation  of  the  camera  optic  axis  in  azimuth  and  elevation.  These 
relations  also  allow  determination  of  the  precise  camera  focal  length 
and  the  rotation  about  the  optic  axis (horizonal  tilt)  needed  for  regis- 
tration of  the  star  field  with  its  image. 

The  next  step  is  to  relate  the  digital  film  plane  positions  of  the 
trail  as  viewed  from  two  photographic  sites.  These  positions  maybe 
coordinates  of  the  center  line  of  the  trail  determined  using  an  automatic 
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scanning  microdensitometer,  or  they  may  be  the  coordinates  of  small 
dots  drawn  onto  the  triangulation  negatives  in  the  regions  of  interest. 
This  latter  method,  while  less  accurate,  must  sometime  be  used  on 
more  turbulent  portions  of  a trail  since  the  automatic  scanner  would 
be  incapable  of  defining  or  following  the  trail  centerline.  The  methods 
of  Reference  13  together  with  camera  orientations  and  locations  are 
then  used  to  calculate  a line  of  sight  vector  for  each  of  the  trail  points 
as  seen  from  the  two  triangulation  sites.  The  best  matching  of  site  1 
points  with  site  2 points  is  accomplished  by  minimizing  the  dihedral 
angular  mismatch  between  the  planes  generated  by  each  point  line  of 
sight  and  the  line  through  the  two  sites.  Finally,  the  positions  of 
matched  points  in  altitude,  latitude,  and  longitude  are  calculated. 


Triangulation  accuracy  at  high  altitudes  is  on  the  order  of  1/4  km, 
a distance  large  compared  to  the  object  size  that  results  from  the  film 
grain  on  the  fast  emulsions  needed  to  record  the  relatively  faint  clouds; 
the  grain  clumping  has  a scale  size  of  the  order  of  20pm  (or,  one 
might  say,  the  impulse  response  of  the  film  is  of  this  order).  K-24 
cameras,  with  their  7-inch  focal  length  lenses,  image  a 25  meter 
(1/40)  spot  at  a (typical)  range  of  200  km  for  a high  altitude  release, 
into  ~ 20pm.  It  is  interesting  to  note  in  this  regard  that  the  resolution 
of  the  grainy  film  sets  a limit  of  about  50  m in  the  triangulation 
accuracy  that  can  be  achieved  with  this  particular  lens.  In  the  lower 
stratospheric  region  and  with  slower  films  (less  grain  clumping)  the 
method  has  produced  positional  accuracies  of  ~10  meters.  Some 
problems  associated  with  the  method  are  detailed  in  Ref  17. 

The  horizontal  wind  field  may  be  calculated  by  analyzing  trail 
motion  for  several  frames  typically  separated  by  intervals  of  10  to 
30  seconds.  Velocity  components  for  each  altitude  are  determined 
by  least  squares  analysis  of  trail  position  versus  time.  At  present, 
the  vertical  wind  component  may  be  calculated  only  for  persistent 
features  which  can  be  uniquely  identified  not  only  from  time  to  time 
but  also  from  at  least  two  sites. 
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V 


These  horizontal  winds  may  be  used  to  study  in  further  detail 
atmospheric  turbulence, heating  rates,  and  diffusion.  One  characteriza- 
tion of  turbulence  which  may  be  derived  from  the  atmospheric  tem- 
perature and  the  wind  field  is  the  Richardson  number  R.  which  equals 

2 1 

(the  Brunt  Vaisala  frequency)  divided  by  (the  vertical  shear  of  the 
horizontal  wind)  . Richardson  number  is  a measure  of  the  local 
stability  or  instability  of  the  atmosphere  and  its  critical  value  is 
usually  taken  as  1/4;  Ri  < 1/4  indicates  instability,  Ri  >1/4  stability. 

Estimates  of  heating  rate  may  be  obtained  as  the  energy  depo- 

2 

sition  rate  per  unit  mass  e = 1/2  i<(9v/9z)  where  v - kinematic  vis- 
cosity and  9v/9z  = the  vertical  shear  of  the  horizontal  wind.  Profiles 
of  horizontal  wind  may  also  be  high  pass  filtered  to  remove  low  frequency 
components  revealing  the  "turbulent"  nature  of  the  velocity  fluctuations  v'. 
Utilizing  the  stability  criterion  Ri  heating  maybe  estimated  for  vertical 
segments  of  the  trail  where  R^  remains 

e = A<  v'3  >/i  (6) 

where  i is  the  segment  length  and  A is  a constant  of  order  unity. 

For  shear  dominated  turbulence,  estimates  of  vertical  diffusivity 
may  be  based  upon  measurement  of  the  horizontal  wind  and  shear. 
Applying  the  one  dimensional  turbulent  energy  equation 


K (9V/9Z)2  - K N2  = e 
m'  n 


(7) 


where  K is  the  momentum  transfer  coefficient,  K,  is  the  heat 
m ^ 

transfer  coefficient,  the  vertical  diffusivity  K may  be  estimated 
when 

R = N2/(9V/9Z)2  <<i  (8) 

and 

K.  * K 
n m 


then 


K(Z)  = 


~ 

(9V/9Zr  m 


(9) 
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A final  estimate  of  eddy  diffusion  for  buoyancy  dominated 
turbulence  may  be  accomplished  using  equations  5 and  6 to  yield 


K = < v'  > i (10) 

where  <v  > is  the  average  turbulent  velocity  over  the  vertical  segment 
of  length  t of  the  velocity  profile.  Restrictions  on  are  the  same 
as  for  equation  6.  Dissipation  and  diffusion  derived  from  equations  6 
and  10  are  probably  underestimated  since  v'  is  likely  to  be  under- 
estimated and  the  i used  is  certainly  smaller  than  the  average  scale 
of  the  energy  containing  eddies. 
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SECTION  IV 


ATMOSPHERIC  MODELING 

Modeling  of  the  concentrations  of  atmospheric  constituents  has 
been  undertaken  over  the  past  few  years  in  an  attempt  to  increase 
knowledge  of  the  processes  which  alter  atmospheric  composition  as  a 
function  of  time  and  altitude  (transport  mechanisms,  solar  heating, 
photoionization,  chemical  recombination).  Computer  modeled  profiles 
of  species  may  be  compared  to  observed  concentrations  to  assess  the 
models'  success. 

These  models  have  become  increasingly  complex  as  they  are 
extended  to  account  for  as  many  known  mechanisms  and  even  very 
minor  ionic  species  as  computer  capability  will  allow.  Serious 
questions  have  been  raised  (Ref  18)  about  the  near  universal  applica- 
tion of  the  assumption  of  diffusive  equilibrium  which  dictates  a 
virtually  static  state  for  individual  species  above  the  turbopause. 
Additionally,  modeled  variation  of  turbopause  height  may  greatly 
effect  species  concentrations  above  the  turbopause  (Ref  19).  It  has 
been  demonstrated  (Ref  20)  that  the  neutral  temperature  profile 
strongly  depends  upon  the  eddy  diffusion  coefficient  (heat  transport 
due  to  eddy  conduction),  and  scale  heights  of  species  are  affected  by 
variation  in  K (by  mixing)  and  by  changes  in  production  and  loss  rates 
(by  changes  in  neutral  temperature).  The  general  use  of  time  invariant 
profiles  for  eddy  diffusion  and  temperature  in  the  solution  of  the  system 
of  mass  and  momentum  conservation  equations  must  receive  further 
attention. 

AFGL  model  codes  for  simple  nonreactive  atmospheres  (Argon- 
nitrogen,  Kryton-nitrogen)  were  modified  to  introduce  transport  due  to 
diurnally  varying  vertical  winds.  Early  attempts  to  execute  this  code 
terminated  when  winds  removed  all  material  from  one  or  more  altitudes. 
Several  factors  leading  to  this  condition  were  identified  and  some  of  the 
problems  corrected. 
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The  first  major  problem  identified  as  an  inherent  instability  in 
the  one-sided  differencing  techniques  used  (briefly  described  in 
Appendix  A).  This  was  remedied  by  instituting  an  offset  in  the  location 
of  action  of  velocities  and  fluxes  - number  densities  and  temperature 
refer  to  the  center  of  an  altitude  cell  whereas  velocities  and  fluxes 
refer  to  tne  upper  boundary  of  the  cell.  Successful  introduction  of 
vertical  winds  must  also  be  done  slowly  - the  wind  increasing  to 
maximum  values  only  after  the  iterative  solution  to  achieve  diurnally 
reproducible  species  concentration  profiles  has  executed  several 
diurnal  cycles. 

A simplified  solution  to  allow  variation  of  mean  mass  and 
temperature  (energy  balance  equation)  was  tested  and  found  inadequate. 
The  theoretical  solution  for  the  full  mass  motion  case  was  derived  but 
not  implemented  because  of  the  considerable  code  modification  that 
would  be  necessary.  Efforts  were  also  made  to  reduce  execution  time 
of  the  code  by  optimizing  certain  calculations  and  by  exploring  how  to 
increase  the  time  step  or  the  altitude  increment. 

The  final  diurnally  reproducible  species  number  density  profiles 
obtained  with  the  inclusion  of  vertical  winds  had  differences  and 
inconsistencies  with  experimental  results  probably  attributable  to 
partial  nature  of  the  code  solutions,  since  some  parameters  known  to 
vary  diurnally  are  still  treated  as  time -independent. 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

Present  techniques  of  spectral  analysis  are  an  effective  means 
of  determining  energy  spectra  representative  of  atmospheric  turbulence 
from  photographic  records  of  chemical  releases.  These  spectra 
allowed  calculations  of  turbulent  atmospheric  heating  rates  and  turbulent 
diffusion  coefficients  which  were  used  as  input  to  computer  codes  that 
model  the  behavior  of  atmospheric  species.  Additional  measurement  of 
heating  and  diffusion  were  derived  from  the  horizontal  winds  obtained 
by  triangulation  methods  from  these  same  photographs.  Large  volumes 
of  data  were  analyzed,  which  were  reported  in  whole  or  in  summary 
form  to  AFGL  scientists  as  the  results  were  obtained. 

Further  improvements  in  modeling  procedures  depend  upon 
incorporating  an  energy  balance  equation  to  allow  use  of  time  dependent 
profiles  for  temperature,  eddy  diffusivity  and  vertical  winds;  extension 
of  the  codes  to  three  spatial  dimensions  must  also  consider  mean 
horizontal  wind  motions.  Adequate  assessment  of  the  performance  of 
more  comprehensive  models  incorporating  these  time  varying  quan- 
tities requires  further  experimental  investigation,  especially  of  eddy 
diffusivity  about  which  the  least  is  known.  The  hundreds  of  existing 
films  of  chemical  releases  taken  over  several  years  can  be  utilized  to 
systematically  study  both  altitude  and  diurnal  variations  of  winds, 
heating  rate,  and  eddy  diffusion  especially  at  altitudes  below  120  km 
where  turbulent  mixing  and  heating  predominate. 
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APPENDIX  A 
MEMORANDUM 

DATE:  June  1975 

TO:  T.  Keneshea 

FROM:  P.  I.  Richards 

SUBJECT:  Instability  of  One-Sided  Differencing 

In  simple  terms,  a set  of  equations  such  as 
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Leads  to  relations  such  as 
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if  one-sided  differences  are  used. 


Now  suppose  a situation  like  the  following  develops 
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where  n.  is  very  much  larger  than  n..  This  leads  quite  properly 
J 4*  1 J — 

i + 1 

to  a large  negative  Uj  , BUT  this  Uj  applies  to  cell  j,  and  its 

l j.  i 

effect  is  to  further  empty  cell  j.  If  Uj  1 in  cell  j - 1 is  not  equally 

large  and  positive  (as  it  may  well  not  be),  n.  decreases  still  further, 
f +2 

making  the  next  U^  even  more  negative,  and  the  situation  runs  away. 
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